Unidirectional d-wave superconducting domains in the two-dimensional t-J model 



O 

o 

(N 
•*-> 
O 

O 



U 

CO 

-4-* 



I 

C 

o 
o 



> 
00 
00 

q 

od 
o 
i> 
o 



Marcin Raczkowski, 1 Manuela Capello, 2 Didier Poilblanc, 2 Raymond Fresard, 3 and Andrzej M. Oles 1,4 

Marian Smoluchowski Institute of Physics, Jagellonian University, Reymonta 4, PL-30059 Krakow, Poland 
2 Laboratoire de Physique Theorique UMR 5152, CNRS and Universite Paul Sabatier, F-31062 Toulouse, France 
3 Laboratoire CRISMAT, UMR 6508 CNRS-ENSICAEN, 6 Bid. du Marechal Juin, F-14050 Caen, France 
* Max- Planck- Institut fur Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 

(Dated: February 1, 2008) 

Motivated by the recently observed pattern of unidirectional domains in high-T c superconductors 
[Y. Kohsaka et at, Science 315, 1380 (2007)], we investigate the emergence of spontaneous modula- 
tions in the d-wave superconducting resonating valence bond phase using the t-J model at x = 1/8 
doping. Half-filled charge domains separated by four lattice spacings are found to form along one 
of the crystal axis leading to modulated superconductivity with out-of-phase d-wave order param- 
eters in neighboring domains. Both renormalized mean-field theory and variational Monte Carlo 
calculations yield that the energies of modulated and uniform phases are very close to each other. 

PACS numbers: 74.72.-h, 74.20.Mn, 74.81.-g, 75.40.Mg 



Puzzling properties of the high-T c superconductors 
have often been attributed to competing instabilities^ 
Indeed, it is believed that doping an antiferromagnetic 
(AF) Mott insulator (which could be described e.g. by 
the so-called t-J model 2 -) leads to quantum disordered 
states with short-ranged magnetic correlations between 
5 = 1/2 spins and exotic properties^ Among them, the 
resonating valence bond (RVB) state was the first theo- 
retical proposal supposed to capture the essence of high- 
Tc superconductivity.— Remarkably, this approach based 
on the Gutzwiller-projected BCS trial wavefunction, the 
parameters of which are usually determined either by us- 
ing renormalized mean-field theory^ (RMFT) or by vari- 
ational Monte Carlo (VMC) method^ not only correctly 
predicted the d-wave symmetry of the superconducting 
(SC) order parameter^ but reproduced in addition the 
experimentally observed doping dependence of a variety 
of physical observables in the SC regime^ 

In fact, exactly at half- filling where a particle-hole 
SU(2) symmetry applies, the RVB phase is equivalent 
to the staggered flux (SF) stated a projected Slater de- 
terminant build from a tight binding model under a stag- 
gered magnetic flux. Remarkably, short-ranged stag- 
gered orbital current correlations have been seen in the 
Gutzwiller-projected <i-wave RVB phased and in the ex- 
act ground state of a small t-J cluster^ Upon dop- 
ing the SU(2) symmetry is broken leading to two dis- 
tinct phases, a d-wave RVB superconductor and a doped 
SF phase, a candidate for the pseudogap phase ; 12 : 13 
characterized by the opening of an antinodal gap in 
the excitation spectrum. Indeed, coexistence of sharp 
nodal quasiparticlcs and broad antinodal excitations have 
been found in angle-resolved photocmission spectroscopy 
(ARPES) studies on an array of underdoped cuprates 
such as La 2 _ :c Sr :!; Cu04 (LSCO)^ Bi 2 Sr 2 CaCu 2 8+5 
(Bi2212)^and Ca 2 _ x Na x Cu0 2 Cl 2 (Na-CCOC)^ 

However, there are also some low-T properties of the 
SC state which cannot be simply explained within the 
original RVB framework. For example, following its the- 
oretical prediction^ static charge and spin stripe or- 



der has been detected in neutron scattering experiments 
and resonant soft :r-ray scattering in some cuprate com- 
pounds as Nd-LSCO±£ and La 2 _ x Ba :c Cu04 (LBCO)^ 
More microscopic evidences of inhomogeneities have been 
given recently by scanning tunneling microscopy (STM) 
probing the doped-hole charge density with atomic 
resolution^ In the SC regime of two different cuprate 
families, Na-CCOC and Dy-Bi2212, bond-centered charge 
patterns with a width of four lattice spacings have been 
seen as well as with two distinct types of spectral gaps: (i) 
a relatively small one on arrays of two neighboring chains 
(or "ladder" ) with clear coherence peaks as expected for 
a d-wave superconductor, and (ii) a broader pseudogap- 
like one on the other separating ladders. This implies 
inhomogeneous superconductivity across the stripe unit 
cell. Moreover, spin-glass behavior has been seen in the 
charge ordered state by muon spin rotation 21 - with a char- 
acteristic length of 2 nm, matching perfectly the STM 
data^ Also puzzling is the linear vanishing of the den- 
sity of states at four nodal points that has been observed 
in the normal phase of LBCO at doping x = l/8»22 

From the above discussion, spatially modulated charge 
order coexisting either with inhomogeneous AF or SC 
correlations appears to be ubiquitous among various fam- 
ilies of underdoped cuprates. While inhomogeneous AF 
order was predicted in the "stripe" scenario^! SC phase 
coexisting with charge order is a new fascinating issue. 
Spatially oscillating d-wave superconductivity appear- 
ing on top of an AF stripe modulation has been pro- 
posed recently^ 3 - We provide here a different theoretical 
proposal based on a RMFT designed to describe quan- 
tum disordered states^ 4 - (no incommensurate AF order is 
then assumed). Although checkerboard states were first 
proposed^ we focus here on modulated SC solutions. 
Our conclusions are supported by VMC data^ 

We start with a t-J model: 



H 



,c ja + h.c. 



-* E (4* 

where the sums run over the bonds. Next, we replace 
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FIG. 1: (color online) Spatial modulation of the hole density 
rihi (labeled from to 7) found at hole doping x = 1/8 and 
t/J = 3 for: (a) bond-, and (b) site- centered 7rDRVB configu- 
rations. Circle diameters are proportional to the hole densities 
while widths of bond lines connecting them are proportional 
to the magnitudes of the pair-order parameter Ay. 

the local constraints that restrict the fermion creation 
operators cj ff to the subspace with no doubly occupied 
sites by statistical Gutzwiller weights^ while decoupling 
in both particle-hole and particle-particle channels yields 
the following mean-field (MF) Hamiltonian,— 12& 

H M F = ~ t E 9ij(4,<r C j,a + h ' C -) -M5Z^,<7 
~ \ J E 9ij(Xji4,<r C j,<r + h - C - ~ \Xij\ 2 ) 

- 1/ E 9ij( A AA-^ + h - c - - i^f). ( 2 ) 

with the Bogoliubov-de Gennes self-consistency condi- 
tions for the bond- Xji = ( c ] a c i a) an d pair-order Aji = 
(cj _ rT c i a ) = (c i _ a Cj cr ) parameters in the unprojected 
state. Moreover, improved Gutzwiller weights depending 
on local hole densities nhi — 1 — X)o-( c l a c i a) are used, 

J = 4 (! -nhi){l -n hj ) 

9ij a l3 + Sn ht n hj (3- 3 (2) + 16f3± (4) ' 

t I An hi n hj {\ - n hi )(l - n h] ) 
9ij V ay + 8(1 - n hi n h j)\xa I 2 + 16|*y | 4 ' U 

with ay = (l- n 2.)(l- n ^.) and/3±(n) = ±\ Xij \ n , 

while the |xy| and |Ay| terms account for the cor- 
relations of the probabilities between nearest-neighbor 
sites^ Using unit cell translation symmetry^ calcula- 
tions (with hole doping x = 1/8 and t/J = 3) were car- 
ried out on a large 256 x 256 cluster at a low temperature 
(3 J = 500 eliminating finite size effects. 

The two modulated RVB states derived from the par- 
ent d-wave RVB phase found in this study are shown 
schematically in Fig. [1] Hereafter we refer to them as ir- 
phase domain RVB (7rDRVB), as they both involve two 
out-of-phase SC domains (see also Ref . [23h . separated by 
horizontal/ vertical bonds with vanishing pairing ampli- 
tudes, named as "domain wall" (DW), where Ay gains 



FIG. 2: (color online) Same as in Fig. [T] but for the 7rDSF 
patterns. Bond lines with arrows indicate the direction of 
charge currents and their widths are proportional to the mag- 
nitudes of the bond-order parameter Xij', +/ — symbols refer 
to positive/negative flux flowing through each plaquette, the 
magnitude of which is represented by the size of +/— symbol. 

a phase shift of n. Specifically, we have found: (i) bond- 
centered 7rDRVB, with a maximum of the hole density 
spread over two- leg ladders depicted in Fig. [Ha), and 
(ii) site-centered 7rDRVB, where the hole rich regions are 
composed of three- leg ladders as shown in Fig. [lib). I n ~ 
tcrcstingly, out-of-phase adjacent SC domains has been 
advocate d 23 i 30 as responsible for some puzzling decou- 
pling of the Cu0 2 planes in stripe-ordered LBCO. 31 

Moreover, as shown in Fig. [21 we also found closely 
related stripes originating from the SF state. An impor- 
tant characteristics of such patterns (named as 7rDSF) is 
again the existence of DWs which (i) act as nodes for the 
staggered current and (ii) introduce into the SF order 
parameter a phase shift of ir. Hence, the 7tDSF phase 
could be regarded as quantum disordered AF stripesJ^ 
In fact, it differs from the pair density wave checkerboard 
solution proposed to interpret the pseudogap stated not 
only in terms of the spatial symmetry but also in that it 
breaks time-reversal symmetry. 

Typical profiles of the site- and bond-centered normal 
7rDSF and superconducting 7rDRVB phases are shown in 
Fig. [3] and the related MF energies are reported in Ta- 
ble [I] compared to VMC estimates^ Here the locations 
of the DWs (oriented along the y-axis) of the SC and SF 
states have been shifted by 2 lattice spacings (exactly as 
in Figs. [T] and ^ for clarity. Although MF results are 
mostly qualitative, as they are exact for particular small 
clusters only^ 3 - they are believed to correctly reflect the 
local spatial structure of the various states, in particu- 
lar the out-of-phase nature of the respective order pa- 
rameters Af£ = g\ i+a A ii+a and modulated flux 
in neighboring domains as shown in Figs. [3{c,d). The 
hole density profiles in the bond-centered (site-centered) 
stripes show the emergence of two (three) inequivalent 
sites in the unit cell. One important feature of these 
DWs is that they are half-filled, i.e. the integrated 
charge density in the direction perpendicular to the 
DWs corresponds to an average of one hole every eight 
sites. Although the large quantitative differences be- 
tween the magnitudes of the charge modulations (as well 
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FIG. 3: (color online) (a,b) Hole density n^i, (c,d) SC order 
parameter Af^ / modulated flux 4?™ indicated by left/right 
arrows, (e,f) spin correlation Sf , and (g,h) bond charge TV, 
found in the 7rDSF (open symbols) and 7rDRVB (filled sym- 
bols) phase. Left (right) panels depict bond- {site- centered) 
pattern; circles (triangles) in panels (c-h) correspond to the 
x (y) direction, respectively. 

as their opposite signs) in the superconducting and nor- 
mal states might be an artefact of the Gutzwiller/mean- 
field approximations^ the generic trend towards charge 
segregation clearly reflects the competition between the 
superexchange energy (Ej) and the kinetic energy (E t ) 
of doped holes. For example, a reduction of the SC or 
flux order parameters (the latter known to frustrate co- 
herent hole motion^) enables a large bond charge hop- 
ping Tf = 2g\ l+y Re{xis+y} along the DWs, as in the 
usual stripe scenario ) 17 ' 29 at the expense of a reduction of 
the AF correlations Sf = -§ g[ l+x (\Xia+x\ 2 + |A M+2; | 2 ) 
along the transverse bonds. Although earlier VMC stud- 
ies^ suggested that the d-wave RVB state has a lesser 
tendency towards phase separation than the SF, our re- 
sults (including the VMC results^ of Table [TJ clearly 
demonstrate that the modulated phases with commen- 
surate period A = 8 (at x — 1/8) are very competitive in 
energy with the homogeneous ones^ Hence, under some 
circumstances, such as long-ranged Coulomb repulsion 25 
or in-plane anisotropy, the d-wave stripe ordered phase 
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FIG. 4: (color online) MF spectrum _E k of: (a) 7rDSF phase 
and (b) 7rDRVB phase along the nodal Y — M path near the 
S point. Dashed line shows of the d-wave RVB phase. 

might correspond to the global minimum as suggested 
by the experimental data for Nd-LSCOf^ LBCO^ and 
Na-CCOC^ all with a tetragonal structural distortion. 
In addition, the reinforcement of the short-ranged AF 
correlations by the finite flux flowing through plaquettes 
of weakly doped regions in the 7rDSF phase helps to op- 
timize Ej. Consequently, such a phase is characterized 
by a better Ej as compared to the uniform SF state (see 
Table [T| . Hence an increased tendency of the latter to- 
wards phase separation should result from increasing J 
or reducing x, in which case the kinetic energy gain is 
relatively unimportant. 

It is instructive to investigate the low-energy part of 
the MF spectrum E^ (see Fig. @| which may reflect in- 
teresting features of the true quasiparticle spectrum. It 
is well known that the uniform SF phase has a cone-like 
dispersion pinned at the S = (tt/2,it/2) point. However, 
doping pushes the node above the chemical potential 
leading to small "pockets" and it has been recently shown 
that a diagonal 7rDSF phase follows the same trend.— In 
contrast, since the vertical 7rDSF state mixes currents 
with different chiralities, one finds two cone-like features 
located symmetrically around the S point. However, only 
one of them crosses the chemical potential which may 
lead to small electronic pockets weakly shifted towards 
the Brillouin zone (BZ) center. They are reminiscent 
of the electron pockets obtained in a recent slave boson 
study of the t — t' — U Hubbard models 7 - Interestingly, 
the parent d-wave RVB phase and its modulated 7rDRVB 
derivate both show a cone-like linear quasiparticle exci- 



TABLE I: MF kinetic energy E t , magnetic energy Ej, and 
free energy F as well as VMC energy Eyyic for a tilted cluster 
of 128 sites of the locally stable phases (all per site): 7rDSF, 
SF, bond-centered 7rDRVB(l), site-centered 7rDRVB(2), and 
d-wave RVB one. ttDSF(I) and ttDSF(2) phases are fully 
degenerate. 



phase 


E t /J 


Ej/J 


F/J 


-Evmc / J 


ttDSF 


-0.8514 


-0.4269 


-1.2783 


-1.3323 


SF 


-0.8622 


-0.4230 


-1.2852 


-1.3389 


ttDRVB(2) 


-0.8736 


-0.4491 


-1.3227 


-1.3359 


ttDRVB(I) 


-0.8719 


-0.4518 


-1.3237 


-1.3359 


d-wave RVB 


-0.8863 


-0.4784 


-1.3647 


-1.3671 
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FIG. 5: (Color) Real part of the SC order parameter A£ c (Q) 
in: (a) d-wave RVB phase with Q = (0, 0) as well as in the 
bond- centered 7rDRVB phase with: (b) Q = ±(tt/4, 0) and (c) 
Q = ±(37r/4, 0) (scaled resp. by factors 4.5 and 8 for clarity). 

tation spectrum pinned at the chemical potential at a 
location in the BZ close to the cone of the vertical 7rDSF 
phase. We also note that the presence of DWs induces the 
formation of a flat band in the vicinity of the antinodal 
direction that could be related to the ARPES observation 
in underdoped cupratesi 14 ' 15 i 16 

Lastly, in Figs. E)Ja,b) we compare the (double) 
Fourier transform of the SC order parameter A^ C (Q) = 
e iQ -^+^ k - Q / 2 M^-Rj)A?P t0 the uniform RVB 

SC gap. At the Cooper pair momenta Q = ±(7r/4, 0) 
both real and imaginary parts (not shown) of A^ (Q) 
retain much of the <i-wave symmetry with a gap close to 
the nodal directions. Moreover, in contrast to Ref. [23|, 
A| C (Q) is also non- vanishing in the 7rDRVB phase for 



Q = ±(37r/4, 0). Of lower amplitude the latter clearly 
reflects the DW structure through its deviation from d- 
wave symmetry, as seen in Fig. EJc). 

To conclude, we note that the relative stability of the 
modulated d-wavc RVB superconductor should be highly 
sensitive to small perturbations of the microscopic Hamil- 
tonian and, hence, can be relevant to some experimen- 
tal cases such as the bond-centered electronic glass with 
unidirectional DWs in Na-CCOC and Dy-Bi2212 seen by 
STM.— In fact, an experimental support to our modu- 
lated bond-centered superconducting state is provided by 
its strong similarities to the experimental pattern indicat- 
ing the spontaneous appearance of two types of parallel 
ladders on nano-domainsi^ 
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